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Abstract. In standard treatments the Eddington luminosity is calculated by assuming that the 
electron-photon cross section is well described by the Thomson cross section which is gray 
(frequency independent). Here we discuss some consequence of the introduction of free-free 
opacity in the Eddington luminosity computation: in particular, due to the dependence of 
bremsstrahlung emission on the square of the gas density, it follows that the associated absorption 
cross section increases linearly with the gas density, so that in high density environments Edding- 
ton luminosity is correspondingly reduced. We present a summary of an ongoing exploration of the 
parameter space of the problem, and we conclude that Eddington luminosity in high density envi- 
ronments can be lowered by a factor of ten or more, making it considerably easier for black holes 
to accelerate and eject ambient gas. 



INTRODUCTION 

The Eddington luminosity plays a fundamental role in our understanding of accretion 
phenomena and, as a consequence, of the AGN physics and evolution. In turn, it is also 
clear that AGNs played a central role in galaxy formation, as testified by several empir- 
ical scaling relations between the mass of the supermassive BHs at the center of stellar 
spheroids (i.e., bulges and elliptical galaxies) and global properties of the spheroids 
themselves. In standard treatments Eddington luminosity is calculated by assuming that 
the electron-photon cross section is well described by the electron scattering Thompson 
cross section. Here we discuss the modifications to this simple picture by taking into ac- 
count bremsstrahlung opacity also: due to the dependence of bremsstrahlung emission 
on the square of the gas density, it follows that the associated absorption cross sec- 
tion increases linearly with the gas density, so that Eddington luminosity is correspond- 
ingly reduced. It is then expected that this opacity source will be important in the range 
spanned by cold accretion solutions (see. Park & Ostriker 1999). These high-density 
conditions were certainly present at the epoch of galaxy formation, thus a quantitative 
estimate of the reduction of Eddington luminosity as a function of ambient gas density 
and temperature is of direct interest when modeling the early stages of galaxy formation 
(see, e.g., Haiman, Ciotti & Ostriker 2003). In the following discussion we often refer to 
formulae provided in Ciotti, Ostriker & Pellegrini (2003, hereafter COP): eq. (1) there 
will be indicated as COPl, and so on. 



THE EFFECT OF FREE-FREE ABSORPTION 



Let Lbh = bMbrc^ be the bolometric luminosity associated with an accretion rate of 
Mbh, where e is the accretion efficiency and c is the speed of light. The momentum 
transfer to the gas electrons per unit time, frequency and volume at the radius r is given 
by (e.g., see Krolik 1999): 
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where A'^ye(v) is the number of photon-electron interactions as given in COPl, ne(r) is 
the electron number density, h is the Planck constant, and at is the total electron cross- 
section corresponding to the various absorption processes (e.g., see Krolik 1999). For 
simplicity we restrict to the optically thin regime: as a consequence, from now on we 
assume Lbh('', v) = Lbh x /bh(v) (see COP). The net acceleration field experienced 
by the accreting gas with density profile p(r) is given by geffir) = — GMbh/^^ + 
Jq{Ap/ At)dv / p{r) where G is the gravitational constant. By imposing geff = one 
obtains 
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when Oi — Oi — SKq'^/Smlc'^ ~ 6.65 x 10~^^ cm^, the classical expression for LEdd 
is recovered. Here we consider two different opacity sources for the emitted photons. 
The first is due to electron scattering as described by the Klein-Nishina cross sec- 
tion <7kn(v) = OidY^^{x), where x = v/Vj and Vj = mgC^/h (see C0P2). The sec- 
ond is due to free-free absorption: its cross cross section Off is obtained from the 
bremsstrahlung emission formula (per unit volume, frequency, over the solid angle, 
Spitzer 1978, hereafter S78) 
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and the Kirchhoff's law jff(v) = KffBv{T), where By{T) = Sithv^c'^ /[exp (hv/kBT) - 
1] and IcB is the Boltzmann constant. Accordingly, 
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where < gff > is the mean free-free Gaunt factor. We assume that a mass fraction 
X of the accreting material is made of fully ionized hydrogen, and the remaining 
fraction 7 = 1 —X by fully ionized helium. Thus, n^ = Xp/mp, nne = (1 —X)p/4mp, 
«e = {^+X)p/2mp, ni = {3 + 5X)p/4mp, and finally J^j^it^iZf = p/mp. In other words, 
at variance with the standard electron scattering case, Off depends (linearly) on the gas 
density. As in COP, we also adopt a (normalized) spectral distribution made by the sum 
of two distinct contributions 
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fx describes a (normalized) double-slope, non thermal distribution of total luminosity 
Lx, with low and high-frequency slopes and ^1 + ^2, respectively, and break frequency 
Vb (COP9). /uv is a normalized black-body distribution of temperature Tuv and total 
luminosity Luv = (COPIO); thus Lbh = ^x+^uv = (1 + <^)^x- Accordingly, we 
study 
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We now compute the four integrals appearing in eq. (6) for some representative values of 
the parameters entering /bh(v). As in COP, we fix ^1 = 0.9 and ^2 = 0.7, and <^ ~ 1 — 2, 
so reproducing the observed Compton temperature in quasar spectra (see, e.g., Krolik 
1999; Sazonov, Ostriker & Sunyaev 2003). We obtain 
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as a very good fit over the range 0.1 < Vb/Vx < 10. The estimate of /J° c^kn/uv^v is 
trivial: in fact, the 99% of Luv is emitted at hv < IO^bTuv hv-r for any reasonable 
value of Tuv- Thus Okn ~ Ot> and from the normalization condition 
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We then obtain 
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where a = T\]\/T, and the dimensionless integral equals 7 + ^(1 + oc), with 7 = 
0.5772... and ^(x) = dlar{x)/dx. Note that if T\jy = T then the integral equals 1; 
moreover, 7+*I'(l + a) ~ for a ^ and ^(1 + a) ~ ln(a) when a ^ 00. Thus 

we have 
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where temperatures are in Kelvin, and we assumed < gff >~ 9.8 (see S78). The integral 
/(T <^«fxdv is affected by a formal infrared divergency if (as in our case) /x ~ v~^i with 
^1 > — 1 : in fact, cTff ~ for v ^ 0. In order to have a proper physical description we 
integrate only at frequencies V > ;tVp, where 
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Vp=jMl^6.35xlO^/ii±^ s-i (11) 
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is the plasma frequency. Under the condition I <^ X ^ k^T /hVp, the opacity contribu- 
tion comes mainly from the Ray leigh- Jeans regime in eq. (4), and the asymptotic value 
of the integral is 

Thus, if Z = 0.75, Vb/Vx = 0.5, = 2 and Tuv = T we obtain 
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For r>10'^ K the free-free opacity contribution to Eddington luminosity is fully domi- 
nated by the low frequency tail in the non-thermal spectral distribution. For example, a 
reduction factor of Eddington luminosity of the order of 10 is obtained for x — 300 and 
;^ ~ 15, when T — 10^ K and T — 10^ K, respectively. 

Of course, the present analysis leaves open several questions, such as the astrophysi- 
cal status of the high-density gas at low temperatures, the characteristic linear scale over 
which the gas becomes optically thick, and the detailed shape of the non-thermal radi- 
ation distribution in AGNs at (very) low frequencies. All these aspects of the problem 
should be proprerly addressed when applying the present results to numerical models. 
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